Self-organization of vortices in type-II superconductors during magnetic relaxation 
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We revise the applicability of the theory of self-organized 
criticality (SOC) to the process of magnetic relaxation in 
type-II superconductors. The driving parameter of self- 
organization of vortices is the energy barrier for flux creep and 
not the current density. The power spectrum of the magnetic 
noise due to vortex avalanches is calculated and is predicted 
to vary with time during relaxation. 



I. INTRODUCTION 

The magnetic response of hard type-II superconduc- 
tors, in particular magnetic flux creep, is a timely is- 
sue in contemporary research (see for reviewmj) . In 
early 60s a very useful model of the critical state 
was developed ta describe magnetic behavior of type-II 
superconductorsoErci One of the distinguishing features 
of this behavior, observed experimentally, is that the den- 
sity of flux lines varies across the whole sample. This 
model of the critical state remains in use, even though a 
significant progress has been made in understanding the 
particular mechanisms Qf|-a-|magnetization and creep in 
type-II superconductorsETH'El. It has also been noted that 
the magnetic flux distribution in type-II superconductors 
is, in many aspects, similar to a sandpile farmed when, 
for example, sand is poured onto a stageQ'D'Q. When a 
steady state is reached the slope of such a pile is analo- 
gous to the critical current density j c of a superconduc- 
tor. Study of the dynamics (i. e. sand avalanches) of such 
strongly-correlated many-particle systems has led to a de- 
velopment of a new concept, called self-organized critical- 
ity (SOC), proposed originally by Bak and co-workersa. 
Tang first analyzed direct application of SOC to type-II 
superconductorsQ. Latep-numerous studies significantly 
elaborated on this topicE£lE£l. 

In practice, especially in high-T c superconductors, per- 
sistent current density j in the experiment is much lower 
than [the critical current density j c due to " giant" flux 
creepB. The concept of SOC, is strictly applied only to 
the critical state j — j c £5~t3 and it describes the sys- 
tem dynamics towards the critical state. Nevertheless, 
it is tempting to analyze magnetic flux creep in type- 



II superconductors during which the system moves out 
of the critical state, in a SOC context, because ther- 
mal activation can trigger vortex avalanches&Eil. How- 
ever, it was found that modifications of the relaxation 
law due to vortex avalanches are minor and can hardly 
be reliably distinguished in the analysis of experimen- 
tal data. Furthermore, flux creep universality has been 
analyticallyj-demonstrated in the elegant paper by Vi- 
nokur et a/.E3. Universality of the spatial distribution of 
the electric field during flux creep has also been found 
by Gurevich and BrandtE£l. The direct application of 
SOC to the problem of magnetic flux creep thus meets 
a number of serious general difficulties. It is clear that 
critical scaling (power laws for vortex-avalanche lifetimes 
and size distributions) observed in the vicinity of the crit- 
ical stateuu must change during later stages of relaxation 
due to a time-dependent (or current-dependent) balance 
of the Lorentz and pinning forces. 

In this paper we propose a new physical picture of 
self-organization in a vortex matter during magnetic flux 
creep in type-II superconductors. In this approach the 
driving parameter is the energy barrier for magnetic flux 
creep rather than the current density. We show that 
notwithstanding its minor influence on the relaxation 
rate, self-organized behavior may be observed by mea- 
suring magnetic noise during flux creep. 



II. BARRIER FOR MAGNETIC FLUX CREEP AS 
THE DRIVING PARAMETER OF 
SELF-ORGANIZATION 

We consider a long superconducting slab infinite in the 
y and z directions and having width 2w in the x direction. 
The magnetic field is directed along the z axis. In this 
geometry, the flux distribution is one-dimensional, i.e., 
B (r, t) = (0, 0, B (x, t)). As a mathematical tool for our 
analyjsisjve use a well known differential equation for flux 
creepM: 

^ = ( Bv o exp (-U (B, T, j) /T)) (1) 

Here B is the magnetic induc- 
tion, v — VQexp(—U(B,T,j)/T) is the mean velocity 
of vortices in the x direction and U(B,T,j) is the ef- 
fective barrier for flux creep. Note that we adopt units 



1 



with ks — 1, thus energy is measured in kelvin. Since in 
our geometry 4irM = J (B — H) dV we get for the mean 
v 

volume magnetization m — M/V from Eq. 



?£ = -A exp (-U(H,TJ)/T) 



(2) 



where A = Hvq/Aitw. 

It is important to emphasize that we do not modify 
the pre-exponent factor Bvq of Eq. [l] or A oLEq. g, 
as suggested by previous works on SOC (see e.g.Q). Such 
modifications result only in logarithmic corrections to the 
effective activation-enetgy, and they may be omitted in 
a flux creep regimet3~tJ. Instead, we concentrate on the 
details of the spatial behavior of flux creep, barrier U (x), 
as analyzed in detail in our previous worktll. In that work 
Eq. [I] was solved numerically and semi-analytically for 
different situations. We emphasize that, in general, the 
barrier for flux creep depends on magnetic field B, and 
persistent current density j (x) is not uniform across the 
sample (see Fig. 1). Thus, j cannot be used as a driving 
parameter for a SOC model. Instead the relevant param- 
eter is U, which stays constant across the sample. Also, 
since experiments on magnetic relaxation are usually car- 
ried out at constant temperature and at high magnetic 
field, we, can assume U (B,T,j) — U (j). Central results 
of RefB are shown in Fig. 1 using a "collective creep" - 
type dependence U (j) = U Q (B/B ) n {(j c /jf - 1), witn 
n = 5 and /i = las an example, .-see also Eq. a below 
(other models are analyzed in RefJl3 as well and produce 
essentially similar results) . Filled squares in Fig. 1 repre- 
sent the distribution of the magnetic induction B (x) / H 
at some late stage of relaxation (so that j < j c ), the 
solid line represents the normalized current density pro- 
file (note that j c is constant across the sample), and open 
circles show the profile of the effective barrier for flux 
creep U (x)/T. All quantities are calculated numerically 
from Eq. pi The important thing to note is that the en- 
ergy barrier U (x) is nearly independent of x, so that its 
maximum variation SUJs of order of T. As also shown 
from general argumentsEZI, such behavior means that the 
fluxon system organizes itself to maintain a uniform dis- 
tribution of the barrier U across the sample. 

The vortex avalanches are introduced in an integral 
way. An avalanche of size s causes a change in the total 
magnetic moment SM = s. This change is equivalent to 
a change of the average current density Sj = SMj = 7s, 
where 7 = 2c/wV. If the barrier for flux creep is U (j), 
then the variation of current Sj leads to a variation of 
the energy barrier 



SU 



dU 



Sj = 7 



dU 
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(3) 



As mentioned above, maximum fluctuation in the en- 
ergy barrier |<5C/| max is of order of T in the creep regime 
(SU « U). Any fluctuation SU larger than T is sup- 
pressed before it arrives to the sample edge due to ex- 
ponential feedback of the local relaxation rate, which is 



proportional to exp(—U/T), (Eq. [l]). This means that 
only fluctuations SU < T can be observed in global mea- 
surements of the sample magnetic moment. Thus, 
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where we denote as s m the maximum possible avalanche, 
which depends on time via dU/dj. It is worth to note 
that Eq. ^ gives the correct dependence of s m on the 
system size and on temperature. It is clear that in a 
finite system the largest possible avalanche must be pro- 
portional to the system volume. Since it is thermally 
activated, it is proportional to temperature T, consistent 
with our derivation. The characteristic time-dependent 
upper cut-off of the avaianche size was experimentally 
observed by Field et alS3 who studied magnetic noise 
spectra at different magnetic field sweep rates, i.e. at 
different time windows of the experiment. 

Our central idea is that in the vicinity of j c the system 
of fluxons, indeed, exhibits self-organized critical behav- 
ior, as initially proposed by TangLI. During flux creep, 
it maintains itself in a self-organized, however not criti- 
cal state in the sense that it cannot be described by the 
critical scaling. The self-organization manifests itself by 
the appearance of almost constant across the sample U. 
Avalanches do not vanish, but there is a constrain on 
the largest possible avalanche, see Eq. ^. Importantly, 
s m depends upon current density and, as we show below, 
decreases with decrease of current (or with increase of 
time), so their relative importance vanishes. 

In order to calculate physically measured quantities let 
us derive the time dependence of s m assuming a very use- 
ful generic form of the barrier for flux creep, introduced 
by GriessercHl. 
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This formula describes all widely-known functional forms 
of U (j) if the exponent a attains both negative and 
positive values. Forj-p = —1 Eq. || describes the 
Anderson-Kim barriercl; for a — —1/2 the barrier for 
plastic creeper is obtained. Positive a describes collec- 
tive creep barrierso. In the limit a ^ this formula 
reproduces exactly logarithmic barrieiO. An activation 
energy written in the form of Eq. || results in an " interpo- 
lation formula" for flux creepcl if the logarithmic solution 
of the creep equation U (j) = Tln(t/to) is appliedl3 (for 
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For a = 0, a power-law decay is obtained: j (t) = 
jc {to/t) n , where n = T/U . 

Using this general form of the current dependence of 
the activation energy barrier, we obtain from Eq. ^ 
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{J) = Wo{j c 
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As we see, the upper limit for the avalanche size decreases 
with the decrease of current density or with the increase 
of time for all a > — 1. For a < — 1 the curvature 
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(9) 



is negative and largest avalanche does not change with 
current, but is limited by its value at criticality) . In 
this case, self-organized criticality describes the system 
dynamics down to very low currents. On the other hand 
the Kim- Anderson barrier must be always relevant when 
j — > jja, thus our model produces a correct transition to 
a self-organized critical state at j = j c . In practice, most 
of the observed cases obey a > — 1 and s m decreases with 
decrease of current density (due to flux creep). 



III. AVALANCHE DISTRIBUTIONS AND THE 
POWER SPECTRUM 

Before starting with calculation of the power spectrum 
of the magnetic flux noise due to flux avalanches, let us 
stress that the time dependence of s m is very weak (loga- 
rithmic, see Eq. ||). This allows us to treat the process of 
the flux creep as quasi-stationary, which means that dur- 
ing a short time, as required for the sampling of the power 
spectrum, current density is assumed to be constant. In 
more sophisticated experimentsEj the external field can 
be swept with the constant rate, which insures that the 
current density does not change, although j < j c . Ac- 
tually, constant sweep rate fixes a certain time window 
of the experiment t/to oc 1/ (dH/dt). Thus, decreasing 
the sweep rate allows the noise spectra to be studied at 
effectively later stages of the relaxation. 

Once an avalanche is triggered by a thermal fluctua- 
tion, its subsequent dynamics is governed only by inter- 
actions between vortices for which motion is not due to 
thermal fluctuations. Thus, we expect same relationship 
between the avalanche lifetime r and its size s as in the 
case of a sandpile: t (t) oc s a (t) and r m (t) oc (t), 
respectively. Using the simplified version of the distribu- 
tion of lifetimes estimated for a superconductor in a creep 
regime from computer simulations by Pan and DoniachEj, 



p (t) oc exp (-r/r m ) , 



(10) 



and assuming that avalanches of size s and lifetime r 
contribute the Lorentzian spectrum, 



L (ui, t) oc 



1 + (ut)" 



(11) 



the total power spectrum of magnetic noise during flux 
creep is 



S{u) 



Using Eq. nfl we find: 



p(r) L(uj,T)dT. 
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(12) 
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(13) 

xv /r]dr] is the expo- 



nential integral. The power spectrum S (to, t) described 
by Eq. |l^ is plotted in Fig. 2 using a solid line. Since 
there an upper cutoff for the avalanche lifetime at r m , 
the lowest frequency which makes sense is 27r/r m . Thus, 
only frequency domain 27r/r m < lu (p > 1) is important. 
In the limit of large p, the spectral density of Eq. [b| has 
a simple asymptote: 



b (UJ) OC 5- 
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(14) 



where j e « 0.577... is Euler's constant. This simplified 
power spectrum is shown in Fig. 2 by a dashed line. For 
p > 10 this approximation is quite reasonable. The usual 
way to analyze the power spectrum is to present it in a 
form S (lu) oc 1/lo v and extract the exponent v simply 
as v = — din (S) /din (ui). In our case the parameter 
p = LOT m is a reduced frequency, so the exponent v can 
be estimated as 



d\n(S) 
d In (p) 



In (p) - 7 e 



(15) 



This result is very important, since it fits quite well the 
experimentally observ ed |. va lues of v which were found to 
vary between 1 and 2E3E3. As seen from Fig. 2, it is 
impossible to distinguish between real \jio v dependence 
and that predicted by Eq. [l3|at large enough frequencies. 
Remarkably, in many experiments the power spectrum 
was found to deviate significantly from the l/u> v behavior 
at lower frequencies, which fits, however, Eq. [T^ . 

Using Eq. [l3| or Eq. |l4| one can find the tempera- 
ture, magnetic field and time dependence of the power 
spectrum substituting p — u)T m — ws^ and using values 
of s m (H, T, t) derived in the previous section. Specif- 
ically, from Eq. || we obtain that any given frequency 
amplitude of a power spectrum increases with time in 
the collective creep regime, but saturates in the case of 
the logarithmic barrier and remains constant in the case 
of the Kim- Anderson barrier. 

In general, we emphasize that the power spectrum of 
the magnetic noise during flux creep depends on time. 
Since parameter p decreases with the increase of time, 
the exponent v becomes closer to 1 during flux creep. At 
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these later stages of relaxation the effect of the avalanches 
is negligible and magnetic noise is mostly determined 
by thermally activated jumps of vortices with the usual 
(non-correlated) l/tu power spectrum. Thus, the man- 
ifestation of the avalanche-driven dynamics during flux 
creep is noise spectra with 1 /lj 11 and decreasing v (t) 
when sampled at different times during relaxation. This 
explains the experimental results obtained by Field et 
al.t3, who measured directly vortex avalanches at differ- 
ent sweep rates. Those found that the exponent v de- 
creased from a relatively large value of 2 at a large sweep 
rate of 20 Gj sec to a smaller value of 1.5 for a sweep rate 
of 1 Gj sec. This is in a good agreement with our model. 

IV. CONCLUSIONS 

In conclusion, self-organization of vortices in hard 
type-II superconductors during magnetic flux creep was 
analyzed. Using results of a numerical solutionis of the 
differential equation for flux creep, it was argued that the 
self-organized criticality describes the system dynamics 
at j = j c . During flux creep, the vortex system remains 
self-organized, but there is no criticality in the sense that 
there are no simple power laws for distributions of the 
avalanche size, lifetime, and for the power spectrum. The 
driving parameter of the self-organized dynamics is the 
energy barrier U (B,j) and not the current density j, 
as proposed by previous work. Using a simple model 
the power spectrum S (uj) of the magnetic noise is pre- 
dicted to depend on time. Namely, fitting S (uj) to a l/u" 
behavior will result in a time-dependent exponent v (t) 
decreasing in the interval between 2 and 1. 
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Figure captions 

Fig.l Results of numerical solution of Eq. [l]for U (j) = 
Uq (B/Bq) (jc/j — 1) at j < j c - Spatial distribution of 
magnetic induction B (x) / H (filled squires); correspond- 
ing profile of the normalized current density (solid line) 
and the corresponding profile of the effective barrier for 
flux creep U (x) /T (open circles). 

Fig. 2 The power spectrum S (ui,t) described by 
Eq. |l3| (solid line) and approximate asymptotic solution 
(dashed line). 
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Fig . 1 Prozorov and Giller 




Fig . 2 Prozorov and Giller 



